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ABSTRACT 

We explore the shape of the galaxy luminosity function (LF) in groups of different mass 
by creating composite LFs over large numbers of groups. Following previous work using 
total group luminosity as the mass indicator, here we split our groups by multiplicity 
and by estimated virial (group halo) mass, and consider red (passive) and blue (star 
forming) galaxies separately. In addition we utilise two different group catalogues 
(2PIGG and Yang et al.) in order to ascertain the impact of the specific grouping 
algorithm and further investigate the environmental effects via variations in the LF 
with position in groups. Our main results are that LFs show a steepening faint end for 
early type galaxies as a function of group mass/ multiplicity, with a much suppressed 
trend (evident only in high mass groups) for late type galaxies. Variations between 
LFs as a function of group mass are robust irrespective of which grouping catalogue 
is used, and broadly speaking what method for determining group 'mass' is used. We 
find in particular that there is a significant deficit of low-mass passive galaxies in low 
multiplicity groups, as seen in high redshift clusters. Further to this, the variation in 
the LF appears to only occur in the central regions of systems, and in fact seems to 
be most strongly dependent on the position in the group relative to the virial radius. 
Finally, distance-rank magnitude relations were considered. Only the Yang groups 
demonstrated any evidence of a correlation between a galaxy's position relative to the 
brightest group member and its luminosity. 2PIGG possessed no such gradient, the 
conclusion being the FOF algorithm suppresses the signal for weak luminosity-position 
trends and the Yang grouping algorithm naturally enhances it. 

Key words: surveys — galaxies: clusters — galaxies: evolution — galaxies: luminosity 
function, mass function 



1 INTRODUCTION 

It is well known that galaxy properties, such as their 
morphology, star formation rate and col our, vary with 
the density of their local environment (|Dresslerl ll98(J ; 
iLewis et all |2002| ; Ikauffmann et al.1 120041 ). Unsurprisingly, 
this is also seen as a variation of t he galaxy luminosity 
function (LF) with environment (e.g. iFerguson et al.l ll99ll; 



Driver et alj Il99ct [Phillipps et all Il998l: IZabludoff et al. 1 



20021 ; iHogg et alj|2003t iPracev et all 



2000; Trcntham et al. 



20051 : (Croton et al.ll2005h 7 These variations are likely due to 



some combination of the suppression of star formation once 
a threshold local de nsity is reached (e.g. iBalogh at al.|[20o3 : 
iTanaka et al.l [2005) and the significant differences in the 
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inter action and/or merger histories be tween environments 
(e.g. lMoore et al.lll99cilBell et al.ll2006l ). 



In a previous paper (jRobotham et al ] |2006f ) we explored 
the variation of galaxy populations across the mass range 
from small groups up to clusters, in terms of the LF param- 
eters for 'c omposite' groups. Using groups from the 2PIGG 
catalogue (|Eke et al.ll2004al ). we showed that the charac- 
teristic magnitude M* brightened and the faint end slope 
a steepened with increasing group luminosity. That is, as- 
suming that the total luminosity of the group members is 
a reasonable proxy for the total mass of the (group or clus- 
ter sized) halo containing them, more massive halos contain 
both brighter galaxies and larger fractions of faint galax- 
ies. This, and the variations seen separately in the blue 
and red galaxy sub-populations, appeared consistent with 
a generic picture of the build up of the giant cluster galax- 
ies through mergers and the amplification of the dwarf end 
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of the LF as 'field' irregulars 'fall in' to larger groups and 
cluster and are transformed, through truncation of star for- 
mation, into dwarf ellipticals. Most of this transformation 
seemed to take place in quite low luminosity groups, i.e. 
quit e early in the hierarchical build-up of the larger s ystems 
(e.g. lHashimoto fc Oemlerll2000l ; IWilman et ai1l2005h . 

In the current paper, we extend this work (and that 
of lEke et al.ll2004b ) by using different LF estimators, sev- 
eral different ways of dividing the groups by 'size' and two 
different group catalogues. In addition we explore the 'en- 
vironment' in more detail by exploring radial LF variations 
within groups/clusters . 

Section 2 describes the input group catalogues, noting 
in particular the differences in grouping algorithms and the 
influence this may have on the subsequent results. Section 3 
discusses the various (group halo) mass proxies and derives 
the LF parameters for each of our composite groups obtained 
from the 2PIGG catalogue. Section 4 presents the equiva- 
lent di scussion for groups in the catalogue of lYang et al.l 
<|2005bh . Section 5 considers the question of the positional 
dependence of galaxy luminosities within groups and Section 
6 summarises our results. 

We assume a standard cosmology with Qm = 0.3, Aa = 
0.7 and a Hubble constant Ho = 100/ikms - . 



2 GALAXY GROUPS AND GROUP 
CATALOGUES 

Since this work is based on the group/cluster environment, 
it will be useful to quantify the defining characteristics of 
these systems. We define the virial radius to be the ra- 
dius within which the average group/cluster mass (dark and 
luminous) density is 200 times the critical density of the 
universe. These radii typically fall in the range 0.3/i _1 Mpc 
< T2oo < 2/i~ 1 Mpc and t he masses within them span 
1O 12 M < M 200 < 10 16 M Q (|Eke et al.ll2004af l. Of course, 
at the largest and most massive end of the scale these sys- 
tems are typically called clusters, but the cut-off between 
the two classes of environment is largely arbitrary. 

As well as affecting the evolution of their constituent 
galaxies, by playing host to the m ajority of galaxy-galaxy 
interactions |Conselice et alj|2005h . groups themselves un- 
dergo a large amount of dynamical evolution. In particular, 
the virial state of the group varies as a function of time; the 
calculated virial mass can vary hugely c ompared to the true 
mass over different dynamic time scales (|Mamo n] (|2007l ) and 
references therein). Significantly, unless the group happens 
to be at the extreme of an expansion after an early collapse, 
the virial mass usually significantly underestimates the true 
mass for groups of all richness classes. In general, groups 
with a smaller crossing time to age ratio will possess virial 
and actual masses that are closer in value. 

Groups are also key in larger scale hierarchical evolu- 
tion as they merge frequently and rapidly in comparison to 
larger structures, and are the building blocks for cluster en- 
vironments. Such a p rocess of group bu ild-up is a generic 
prediction of ACDM (|Lacev et al.l jl993|) et seq.), and it is 
found that merger rates are roughly proportional to the mass 
of the system, with the consequence that present day groups 
can trace their half-ma ss progenitor groups to a higher red- 
shift than can clusters (Ivan den B osch 2003). 



2.1 Approaches to Grouping Galaxies 

The main trade-off that in general has to be made when 
choosing a grouping algorithm is whether it is more desir- 
able to find all the bound members of a group and suffer 
some interlopers as a consequence, or to minimise the inter- 
loper rate and not group all bound objects. Tied into this 
compromise is the fact there is no correct answer when de- 
termining the strict boundaries of a group; depending on the 
over-density limit applied, the final product will have differ- 
ing degrees of the infall region included in the final system. 
Ma ny different approa ches to grouping galaxy data are used 
(see lGerke et al.ll2005l . for a recent discussion), but the two 
most relevant ones for our work are discussed below. 



2.1.1 Friend- Of- Friend (FOF) 

All forms of the FOF algorithm rely in some way on tak- 
ing a galaxy in the sample and searching the surrounding 
area/volume for other galaxies that meet certain criteria 
for separation in both velocity and projected coordinate 
space. The p i oneer ing application of this process was by 
iHuchra et all (|l982h . where the catalogue of lde Vaucouleura 
( 19751 ). the earliest reasonably complete attempt at a group 
catalogue, was reconstructed using fully quantitative means. 
The 2PIGG catalogue used for this work is an adapted 
FoF algorithm which requires that linking criteria are met 
in both projection and radial separation. A brief overview 
of the 2PIGG algo rithm is given bel ow, but the interested 
reader is referred to lEke et al.l (|2004al ) for a detailed descrip- 
tion of the process. 

When linking is considered, it is always in respect to 
how numerous objects are expected to be at a given redshift 



n{z,6) = 



( p[M{b j ,z)]c z (b 3 ,e)db J 



(1) 



wher e 4>\M(bj,z)] re fers to the Schechter luminosity func- 
tion flSchechterlll976T ) and has been determined for 2dFGRS 
in multiple pap ers; the function used for 2PIGG is that of 
Norb erg (2002). c z is a completeness correction determined 
from the position of a galaxy on a pixelated completeness 
mask, and typically scales the local density by a factor ~ 0.8. 
The integral is from the bright limit of the galaxy sample 



(effectively b 



bright 



= — 00) and the faint limit for each galaxy 



19.4 is the target magnitude for 2dFGRS, but the 



exact depth varies). M(bj , z) uses the redshift of each galaxy, 
and hence luminosity distance, to convert the apparent mag- 
nitude into the absolute magnitude required for the LF. The 
next step is to calculate linking lengths for each galaxy that 
represent a certain degree of over-density relative to n(z, 9): 



h 



b 

777/3 



Rli 



(2) 



where b is the scaling factor required for each galaxy link in 
order to achieve a certain degree of overdensity. Typically, 
N-body simulations opt for a relative over-density 125, so 
that the scaled linking length is b=0.2 times the mean in- 
terparticle separation (since the overdensity is oc 6 -3 ). This 
does not imply the mean galaxy separation factor b ga i = 0.2 
precisely, since dark matter and luminous galaxies (whilst 
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tracing each other well in general) are not identicall y dis- 
tributed. In fact it is predicted by ICole et all ( 2000) that 
light is more concentrated than matter in a galaxy group 
halo, so bg a i is likely to be smaller than 0.2 to achieve the 
same level of dark matter clustering. R is the radial scaling 
desired, and must be larger than 1 in order to account for 
peculiar velocities within group haloes. Thus two galaxies 
i and j, at a comoving distance d c ,i and d c j with angular 
separation 9ij, are linked if 

and 

\d c ,i-d c>j \ < + (4) 

Neither b nor R are fixed, instead they are allowed to vary 
to account for the finding that small groups tend to have 
overestimated sizes and velocity dispersions and vice-versa 
for large groups (i.e. b ga i has a slight dependence on group 
mass). Thus 

6 =^(^r (5) 

and 

* = ^(^£f (6) 

where A is the observed overdensity of a cylinder with a co- 
moving projected radius of 1.5/i _1 Mpc and an aspect ratio 
of Rgai- Applications to mock catalogues determine appro- 
priate values of the parameters to be A fa — 5, tb = 0.04 
and £r = 0.16, and most s ignificantly b pa i = 0.13, R ga i = 11 
and Lxmax = 2/i~ 1 Mpc (|Eke et all 1200481 ). An extra lim- 
itation was applied to the FOF so that no perpendicular 
linking length could exceed 2/i -1 Mpc, and by implication 
22/i -1 Mpc in radial coordinates, though in reality such an 
extreme length will never be the limiting factor for the radial 
distance. 

The combination of these parameters, and the use of a 
cylindrical rather than spherically optimised linking volume, 
was found to have the desired effect of maximising the truly 
grouped objects associated with a dark matter halo (group) 
and, depending on the redshift r anges used, minimising the 
interloper rate (|Eke et al.ll2004ah . This latter factor was the 
main influence in determining a reasonable redshift limit for 
our subset of the group catalogue, as will be discussed in 
more detail later. 

The main issue with using a FOF method is that the 
overdensities are by definition local, in fact since the over- 
densities are based on the nearest neighbour they could 
not be calculated with less data. Background tests carried 
out using Monte Carlo simulations of galaxy distributions 
<|Robothaml 120081) demonstrate that in a volume with ex- 
actly the desired grouping density, the multiplicity of the 
system has a stochastic effect on how effectively the FOF 
algorithm groups all the objects. So even though the aver- 
age density in the volume is precisely what we require, in- 
dividual galaxies can be missed due to positional noise and 
nearby (sub)systems can avoid association. 

After running a large number of realisations (100,000 di- 
vided by the system multiplicity) , it was clear that the lower 



multiplicity systems suffered the former problem most: sys- 
tems with multiplicities 5, 10, 20, 50, 100, 200 and 1,000 
achieved 84%, 89%, 91%, 94%, 96% and 97% grouping com- 
pleteness when using 1.05 o as the FOF linking length, where 
So is the mean interparticle distance within the overdensity. 
As for the latter problem of subgroups not associating, the 
linking length that, on average, will link all objects via a 
FOF routine (essentially the overall percolation length) ac- 
tually increases slightly as a function of multiplicity, chang- 
ing from 1.1 to 1.3S over the same range of multiplicities. 
Thus it is clearly important to consider which volume aver- 
aged overdensity is of interest when constructing a grouping 
algorithm, and whether it is worth increasing the interloper 
rate to ensure all truly associated objects are grouped; in 
the case of 2PIGG, the creators decided the answer to this 
question was yes. Thus FOF parameters were chosen so as 
to ensure nearly 100% of true group members were success- 
fully grouped, with the unavoidable side effect of a non-zero 
interloper rate. 



2.1.2 Halo Based Grouping and The Conditional 
Luminosity Function (CLF) 

Halo based methods for grouping centre on the premise that 
when creating a group catalogue the most important crite- 
rion to achieve is that all the grouped galaxies share a com- 
mon dark m atter halo, the early ve rsion of such a technique 
being that of IPostman et afl (|l996l ) . The main principle be- 
hind such methods is that dark matter halos have a certain 
probabilistic distribution, and then depending on the mass 
of a particular halo the LF populating it is also probabilis- 
tic (th e Conditional Luminosity Function, CLF from here 
on; see Ivan den Bosch et alj|2003l ). Thus it seems intuitively 
clear that we can create a hypothetical Universe going one 
way - building the halos then populating them with galaxies 
- but what is required from a redshift survey is to take the 
galaxy data and try to go the other way and match them up 
with the expected distri bution of da r k matt er halos. This 
is the approach used by lYang et al l l|2005bl ) for grouping 
galaxies in the 2dFGRS and in the Sloan Digital Sky Sur- 
vey (SDSS). 

Their first step was to make a rough FOF catalogue of 
the survey in question using quite restrictive conditions: in 
the notation used in the previous section b = 0.05 (rather 
than ~ 0.13) and R = 6 (rather than ~ 12). The reason 
for such extreme overdensities was to find the likely group 
centres, not complete groups, and Yang et al. demonstrated 
this to be an effective means of locating the likely sites for 
halo centres. On top of the FOF centres, individual galaxies 
were also added to the list of potential centres if they were 
the brightest galaxies within a cylinder of radius l/t -1 Mpc 
and velocity depth ±500 kms -1 . 

Once these potential centres were identified an iterative 
process of assigning galaxies to the halos was undertaken, 
and for a given halo-galaxy assignment calculation several 
variables were considered. The first was the current mass of 
the halo, as determined by the luminosity of all the galaxies 
grouped so far (taking into account luminosity dependent 
incompleteness) L group . This was then extrapolated to find 
the total halo mass using: 
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= L 



group r co 



J °° L<p(L)dL 



J™ L(f>(L)dL' 



(7) 



where <f> is the LF (as above) and Ln m is the minimum 
observable galaxy luminosity at the redshift of the group. 
The same LF function parameters were used throughout for 
this algorithm, but to return the expected CLF, Yang et al. 
used an empirically motivated luminosity dependent mass- 
to-light ratio. Interestingly though, they found that using 
a constant M/L = 4Q0h(M / L)q returned almost identical 
groups, lending credence to the robust nature of their ap- 
proach. 

Based on these masses for the group halos, the group 
velocity dispersions and vir ial sizes could be estimated us- 
ing the CDM halo density of lNavarro. Frenk fc White! i| 19961. 
NFW hereafter) 



p(r) 



Po 



(r/r s )(l + r/r s ) 2 



(8) 



where r 3 is the characteristic radius and p is the density 
when r/r a ~ 0.47. Once dispersions and sizes had been es- 
timated, any potential addition to the group was checked 
against the criterion that it was in a region of overde nse 
phase-space relative to the group (see lYang et alJfeoOSbl for 
details) . 

When all possible additions to all of the groups had 
been made, the group properties were recalculated and the 
whole process repeated to see if the new groups could add 
any members that previously failed to meet the phase-space 
overdensity criteria. This algorithm was applied iteratively 
until the catalogue was stable and no new galaxies were 
grouped. 

The major difference between the purely FOF approach 
to grouping and the slight mixture of techniques used for 
the halo based grouping is the scale at which a potential 
group member is considered. In the FOF algorithm purely 
galaxy-galaxy relations are calculated and by inference from 
this, groups of associated objects are constructed, leading to 
the possibility of incompleteness as discussed earlier. In the 
case of halo based grouping, the entire system built so far 
is used to ascertain grouping of any further galaxies; the 
only variables are the relative projected distance and radial 
velocity displacement from the currently constructed group. 
In comparison, the FOF approach does not vary the linking 
length as a function of what has already been grouped, thus 
a massive overdensity on one side of a galaxy does not extend 
the linking length, even though it intuitively should. The 
qualitative advantage of such an approach is that it requires 
less prior information in order to construct a catalogue; both 
approaches to group finding require knowledge of the LF, 
but only the halo based method requires knowledge of the 
density profile — empirically well understood for clusters, 
but less so for small groups. 



2.2 Comparison of the 2PIGG and Yang-2dF 
Catalogues 

As discussed in detail above, the 2PIGG catalogue of 
lEke et al l (|2004al ) uses a FOF algorithm, calibrated against 
a test mock galaxy survey, on galaxy data provided by 
the 2dFGRS. Because of restrictions imposed by Eke et 
al. only 191,440 of the final 245,591 2dFGRS galaxies were 



used, though this is significantly more than in the other 
three 2dFGRS based group catalogues of Yang et alj||2 005b1 
henceforth the Yang-2dF catalogue: iMerchan et all ( 2002) 
which is not publicly available and iTago et al.l (|2006T) . Of 
the known cluste rs in the volume concerned in the Abell 
llAbell et al.lll989h. EDCC ijCollins et alJl995l ) and APMCC 
l|Dalton et alJll99Tt > samples, the 2PIGG algorithm returns 
51% (84 of 166 clusters with a known redshift in SIMBAD) 
where the match distance in comoving space between iden- 
tified group/ cluster centres must be within 2/i _1 Mpc. The 
same match requirements for Yang groups only returns 28% 
(47 of 166). 

Compared to the other catalogues, 2PIGG has a large 
number of groups where the number of members is > 4; 
7,020 groups (comprised of 55,753 from 191,440 possible 
galaxies) meet this criterion, compared to 2,209 groups 
(from 60,000 galaxies) in Merchan & Zandivarez, 2,471 
groups (from 150,715 galaxies) in Yang-2dF, and 3,044 
groups (from 184,495 galaxies) in Tago et al. Proportion- 
ately 2PIGG is closest to the Merchan & Zandivarez cata- 
logue for this figure. In terms of galaxies being grouped (2 or 
more members) or not, 2PIGG is most similar to Yang-2dF, 
where 60% are grouped compared to 2PIGG's 55%. Both 
these figures are somewhat different to the 30% obtained by 
Tago et al.; Merchan & Zandivarez only considered groups 
of 4 or larger so cannot be used for this last comparison. 

A major difference between Yang-2dF (and also the 
Tago catalogue) compared to 2PIGG is their respective 
treatments of group substructure; in particular the treat- 
ment of single galaxies and halos. A single galaxy can be 
considered as belonging to a small group, or alternatively 
it could be considered to be part of the group halo as rep- 
resented by just o ne galaxy. This prob lem of substructure 
was considered bv lBurgett et all l)2004h . the conclusion be- 
ing that most groups have substructure that will merge over 
long dynamic times. The problem faced when grouping is 
whether to count these merging groups as one group or 
as many sub-groups. Yang et al. biased their algorithm to- 
wards more compact, fully formed and gravitationally bound 
groups, hence the significantly larger fraction of single galax- 
ies in their results. Both methods can be justified, but it is 
worth considering why the results appear to be so different. 
To an extent it is a matter of taste when considering this 
problem, but we should be aware later that the results we 
obtain using the 2PIGG catalogue are based on groups which 
may be in the merging process rather than fully formed (at 
least relative to those in the Yang-2dF sample). The main ef- 
fect of this will be realised for low multiplicity groups, where 
the amount of clustering may seem overemphasised. 

Figure[T]contains two plots that demonstrate significant 
differences between the two group catalogues used in this 
work. For the full volume limited samples it is clear from the 
left panel that 2PIGGs have a higher typical velocity disper- 
sion. The right panel compares velocity dispersions between 
groups that best match each other. The best match is de- 
fined as the two-way match that has the largest product of 
the galaxy matching fraction for each group. An example of 
matching groups would be a case where 4 galaxies are found 
in both a Yang-2dF group and a 2PIGG and the Yang-2dF 
group contains in total 5 galaxies and the 2PIGG contains 7. 



The product of these fraction is 
native matches might be, say, i- 



0.46. 



0.2 and §§ 



The best alter- 
0.43, so 
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the match described first would be the preferential match. 
When reference is made to matched groups, it is in the sense 
described here. Obviously it is possible to devise alterna- 
tive methods of group matching, e.g. based on spatial close- 
ness or brightest galaxy, but for illustrative purposes the 
approach described is both simple to calculate and inter- 
pret. The left panel of Figure [2] demonstrates separately for 
Yang-2dF groups and 2PIGGs what fraction of each group 
is contained within its best matching counterpart. It is in- 
teresting to note that a much larger number of Yang-2dF 
groups are 100% contained with the best matching 2PIGG. 
Also, a significant fraction of 2PIGGs share less than half 
their galaxies with the best matching Yang-2dF group, again 
indicating that on occasion multiple Yang groups are merged 
to form a single 2PIGG. 

Based on this matching, the right panel of Figure [T] 
displays an unambiguous tendency for 2PIGGs to possess 
much larger velocity dispersions for groups that are trying 
to describe the same halo. The median velocity dispersion 
of groups with multiplicity > 5 (calculated using the same 
GAPPER algorithm described in Eke et al. 2004a) is 199 
kms" 1 for Yang-2dF groups and 238 kms -1 for 2PIGG. In 
comparison, the best matching groups have median veloc- 
ity dispersions of 206 kms -1 for Yang-2dF groups and 282 
kms for 2PIGG. This suggests that the groups created by 
2PIGG that have no counterpart in the Yang-2dF catalogue 
are of lower velocity dispersion, since the median velocity 
dispersion is so much higher for matched 2PIGGs. 

The most significant difference, in terms of how it might 
affect the present work, is the general tendency for the 
2PIGG algorithm to link more objects compared to the Yang 
et al. group finding algorithm. Using the same redshift limits 
of 0.05 < z < 0.1 and a minimum of five galaxies in a group 
(see below), 2PIGG assigns 17,805 galaxies to 1,535 groups 
(average of 11.6 per group) compared to 5,364 galaxies in 
544 groups (average of 9.9 per group) in Yang-2dF. Over 
the full multiplicity range, the Yang-2dF catalogue has a 
slightly higher proportion of low multiplicity groups (though 
less in total), and also has an upper multiplicity limit of 89 
compared to 2PIGG's largest group of 163. 

Inside the range considered in this work, 77 of the 
halos are found to be identical (comprised entirely of the 
same galaxies). In total 327 groups are unambiguous bijec- 
tive matches, meaning more than 50% of galaxy members 
in both the Yang group and 2PIGG can be found in the 
best two-way matching halo. Considering the Yang groups 
alone, 532 of the 544 groups considered have a 2PIGG sys- 
tem that contains more than 50% of their assigned galaxies. 
This large fraction is possible because 47 2PIGGs contain 
multiple Yang groups that meet this criteria. The largest 
2PIGG, a group comprised of 165 galaxies, contains 6 Yang 
groups that owe at least 50% of their members to this single 
large 2PIGG. 

In total we find 3,572 galaxies in the Yang-2dF groups 
from this sample of similar groups, compared to 6,383 galax- 
ies in the 2PIGG groups. The right panel of Figure [2] com- 
pares the multiplicities between these similar groups found 
in the two different catalogues. There is a clear tendency for 
the 2PIGG catalogue to group additional galaxies compared 
to the Yang-2dF catalogue, in fact only 10 similar groups are 
larger in Yang-2dF. By the large quantity of bijective match 
halos we can be confident that they are, in general, describ- 



ing the same dark matter halo, the main variance being the 
degree to which galaxies are found to be associated with it. 
We can conclude that the 2PIGG catalogue describes the 
same base population for a given halo as the Yang-2dF cat- 
alogue, but it also tends to include a more tenuously asso- 
ciated population, which may be less bound to the halo in 
question. 



COMPOSITE 2PIGG LUMINOSITY 
FUNCTIONS 



As in put for our calculations we follow iRobotham et all 
i|2006h in first selecting all groups from the 2PIGG cata- 
logue (|Eke et al.1 [2004a! ) which have 0.05 < z < 0.10. The 
z — 0.05 lower limit is motivated by the small volume sam- 
pled by 2PIGG at low redshift (so that groups may not be 
representative), while at z = 0.10 the 2dFGRS apparent 
magnitude limit begins to exclude even moderate luminos- 
ity galaxies (Mb < — 18 + 51og/i) from the sample. [Note 
that hereafter the term 5 log h will be suppressed for clarity, 
so unless otherwise indicated, absolute values will be quoted 
as if h = 1.] In addition, over this ra nge of redshift th e in- 
terloper rates remain low and steady l|Eke et al.ll2004aT ). 

Since 2PIGG includes groupings as small as 2 galax- 
ies, we also impose a minimum group size of 5 members. 
(For comparison, at the range of distances we use, the Lo- 
cal Group would appear to have a multiplicity of 2 to 4). 
This should increase the likelihood that the groups we se- 
lect are actual bound physical systems. As discussed in 
IRobotham et al.l(|2006h . these choices leave 1535 groups con- 
taining a total of 21,752 galaxies, when corrected for incom- 
pleteness. The incompleteness arises primarily as a result 
of fibre collisions in the 2dFGRS , that is two fib res have 
a minimum possible separation. lEke et ail (|2004bl ) present 
details of the calculation and we follow their prescription, 
weighting each observed galaxy by a factor depending on 
its local galaxy density (as determined from the positions of 
its 10 nearest neighbours). These groups can then be used 
in various combinations to determine composite luminosity 
functions for galaxies in different environments, as described 
in the following subsections. 

To derive the LFs, we first calculate absolute Mf, . mag- 
nitudes and rest-frame bj — tf colo urs for galaxies in the 
groups, following the prescriptions of ICole et al.l (|2005l ). We 
ignore the relatively small (and uncertain) amount of evolu- 
tion within the narrow redshift range we survey. In addition 
to the LFs derived for the entire sample we also create LFs 
for red (quiescent) and blue (star-forming) ga laxies sepa- 
rately. For this we again follow lCole et all <|2005l ) and divide 
the samples at a colour of bj — tf — 1.07. 

Before attempting to determine any observed trends, we 
should consider the way we quantify the group or cluster en- 
vironment. A group's total mass is evidently a critical mea- 
sure, as halo mass is the key variable in structure formation 
models (Press & Schechter 1974 et seq.; Y ang et al. 2005c). 
The t otal blue luminosity, M^. , as used in IRobotham et all 
(2006) has been demonstrated to be a good measure of (and 
therefore proxy for) total stellar mass and total halo mass 
(Padilla et al. 2004, Yang et al. 2005d, 2005e), but one 
could argue that the total R band luminosity, M^ F , would 
be preferable to the B band luminosity, as it is less sub- 
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Figure 1. Plots comparing the velocity dispersions for 2PIGG and Yang-2dF groups. The left panel is a comparison PDF for all systems 
with multiplicity larger than 4 and 0.05 < z < 0.1. The solid vertical line indicates the median for each distribution, and the dotted lines 
are the inter-quartile range. The right panel shows the different values for the velocity dispersion between best matching groups. It is 
quite evident that in general 2PIGGs have greater velocity dispersions and a broader distribution of velocities. This is the case evident 
for matching systems also, so it is not merely an effect of the 2PIGG algorithm creating many large velocity dispersion groups that the 
Yang algorithm misses. 




Figure 2. Left panel describes what fraction of each 2PIGG/ Yang group can be recovered in its best matching (in terms of two-way 
membership) Yang/ 2PIGG counterpart. A large number of Yang groups are 100% contained within a 2PIGG, but the opposite is 
generally not the case. Right panel is a comparison of group multiplicities for best matching group halos. The background grey-scale 
image is a guide to the local log-space density of objects, where darker regions reflect a larger number of groups. 



ject to contributions from short term star formation. In fact 
Robotham (2008) found essentially the same results when 
ordering groups by red luminosity as Robothem et al. (2006) 
did when using total bl ue luminosity. We refer the reader to 
iRobotham et all l|2006l ) for details of our original results and 
for a comparison of our and earlier results on the variation 



of the LF with group luminosity (e.g. Eke et al. 2004b, Zan- 
divarez, Martinez & Merchan 2006). 

Because different environmental indicators may each 
have their own merits, in the following, we use two different 
methods for sorting the data, a virial mass estimate, and the 
multiplicity (number of members). In general,these parame- 
ters all correlate with each other, and with group luminosity 
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Figure 3. Comparison of different methods for determining the relative size of galaxy groups. Red lines show the best fit linear trends 
in the data as a visual guide. The trend shown in right panel is fixed such that five magnitudes spans two decades in mass — as would be 
expected if the relationship between mass and luminosity scaled linearly. 



(see figure |3}, but at a finer level they show disagreements 
as to which groups should be considered most similar. 



3.1 Method 

In lRobotham et all (|2006l ) we computed a composite ga laxy 
LF, following the method described by ICollessI (|l989h . To 
ensure the veracity of our results, for the present paper, we 
instead use a maximum likelihood (ML) method to fit the 
LFs for the composite groups. 

A key improvement when using the ML method is that 
the data is not binned, instead a probability density func- 
tion is constructed based on the Schechter luminosity func- 
tion itself. To do this correctly the integrated luminosity 
function is scaled to 1 for each observation, so for any given 
galaxy this scaling factor is determined by the integral be- 
tween — oo and the absolute magnitude limit of the survey 
for that galaxy — this is easily found by taking the appar- 
ent magnitude limit of the survey and subtracting the k- 
corrected distance modulus for the galaxy. Once the correct 
scale has been calculated, the log of the density of the PDF 
at the magnitude of the galaxy is determined and then this 
whole process is repeated for every galaxy in the subset of 
interest. The sum of the logs is finally calculated, giving 
us the 'log-likelihood', and then the Schechter parameters a 
and M* are repeatedly changed. The combination of param- 
eters that best fit the data, or strictly speaking the fit that 
yields the 'most likely' data, is that which maximises the 
log-likelihood. The confidence limits in log-likelihood space 
are much like those of x 2 space except compacted by a fac- 
tor of 2, i.e. for the error about one parameter (so leaving 1 
degree of freedom) the 95% confidence limit is found at the 
log-likelihood value 1.92 less than the maximum, instead of 
3.84 as would be expected when integrating the \ 2 distribu- 
tion. 

An advantage of using the maximum likelihood tech- 
nique here is that it does not require the completeness at 
the faint end to be corrected for on a group by group basis, 
instead it is taken into account for every galaxy by virtue 
of the magnitude limit of the survey which is used to scale 
the PDF. This is somewhat akin to a 1/Vmax approach to 
fitting, but due to the galaxies not being binned, extreme 
distortions at the faint end are not likely to occur. Thus 



this approach, whilst agreeing excellently with the previous 
paper in like-for-like systems, is not prone to binning errors. 

To determine the quality of fit to the data, if we as- 
sume the correct PDF has been applied to the data, then 
it would be expected that the distribution of the integrals 
of the Schechter luminosity function between — oo and Mt >j 
for each galaxy should be uniformly distributed. The degree 
of uniformity for this distribution can be determined in a 
number of ways, but here the KS-test has been used. The 
p-value returned from the test will give us an indication of 
the quality of fit. This means that as well as obtaining the 
best fit given an assumed Schechter model, we also have a 
means of quantifying whether the model is appropriate. 



3.2 Virial Mass 

For the mass ordered data, we follow Robotham et al. (2006) 
and again use 10 bins with as near to equal numbers of 
galaxies as possible. Figure [4] contains the variations of a 
and M* as a function of virial mass as determined by the 
method of Eke et al. (2004b), using the radius (specifically 
the root mean square radial distance of galaxies from the 
group centre, r rms ) and velocity dispersion parameters in- 
cluded with the group catalogues. This assumes a degree of 
excess KE compared to GPE in the relative virialisation 
(-2KE/GPE = 1.2), but was found to correspond well to 
the underlying group mass. 

Both the total population and the red galaxy popula- 
tion show evidence of M* becoming brighter as a function 
of group mass, supporting the previous results, while the 
blue galaxy trend appears to be more complicated, with M* 
first becoming fainter with group mass, and then brighten- 
ing sharply for the most massive groups. However, for the 
sub-populations there is more scatter than with the lumi- 
nosity orderings, so, for instance, M* originally becoming 
fainter for the blue galaxies is quite possibly the result of 
an outlier at group mass ~ 10 14 M Q (given the suspiciously 
large a value at this point and the well known correlation of 
errors in M* and a). The larger scatter is likely due, at least 
in part, to the fact that the calculated virial masses are not a 
particularly good indica tor of mass until we reach the m ost 
massive groups (see also lEke et al 1 l2004bl : lMam"ol]|2007l ). Of 
course, in the virial mass case there is no (direct) selection 
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against bright galaxies being in small, low mass groups (as 
there is for low luminosity groups), so any such cases might 
also increase the scatter. On the other hand the continued 
existence of the same trend, particularly for the red galax- 
ies, confirms that these were not forced by the luminosity 
selection in the previous work. 

The trend seen in a values is more straightforward, all 
galaxies combined show a subtle but recognisable pattern 
of q becoming steeper as a function of group mass, start- 
ing and finishing with similar values to the extremes seen in 
Robotham et al. (2006) and Robotham (2008): -0.9 for the 
least massive up to —1.2 for the most massive. When this 
is subdivided, the red galaxies show a very strong steep- 
ening trend as the faint end becomes systematically more 
populated as a function of group mass, whereas the blue 
gala xies are consistent w i th an almost constan t a ~ —1.2 
(c.f. iPopesso et all 120061 ; iMercurio et~afl l2006h . with only 
the most massive systems displaying a notably steeper 
faint end (a ~ —1.5) , as frequent l y rep orted for clusters 
(e.g. iDriver et al.lll994T ). iBell et al] (2004 ) and ICowie et all 
1 2008) note a similar large change in the shape of the mass 
function of red galaxies in general, compared to that of blue 
galaxies, due to the steady migration of 'blue cloud' galaxies 
to the 'red sequence'. The early type LF thus tends to be 
affected at the fain t end due to the recent creation of dEs 
|Smith et alj|2008h . while once the late types ( of all lumi- 
nosities) feel the affect of the dense environment they quickly 
cross the 'green valley' and are removed from the late type 
sample measured here. This is, of course, consistent with 
the recent build-up of the red sequence after galaxies hav e 
fallen into groups and clusters (see, e.g.. | Stott et al.l (|2007h ; 
De Lucia etal] d2007f); ISmith et all d2009h for clusters, and 

"l2ooj ; r" 



Tanaka et al.l (|2005l ); iHaines et all 120081 ) for lower density 



environments) . 

The observed behaviour of the mass ordered groups con- 
firms and strengthens the agree ment with the pred ictions 
of the halo occupation model of I Yang et al.l (|2005al ) noted 
by Robotham et al. (2006). At low group mass (luminos- 
ity) Yang et al. find that the luminosity of the brightest 
galaxies (so, effectively, M*) should increase quite steeply, 
with a much slower increase in the more massive groups, 
the turnover in sl ope occuring at halo masses a few times 
10 13 M n (see also [Vale et all |2004| ; ICoorav fc Milosavlievicl 
120051 ). lYang et alj ""(f2005a) argue that this change in be- 



haviour occurs at the transition from efficient to inefficient 
cooling in different mass halos, and is the same feature re- 
quired in semi-analytical modelling of galaxy formation in 
order to match the bright end of the overall galaxy LF (e.g. 
iBenson et all 120031 ). while Cooray & Milosavljevic suggest 
that it marks the mass at which the dynamical friction tim- 
scale for orbital decay in the halo matches its age. Observa- 
tionally, Eke et al. (2005) find that the variation of group 
mass-to-light ratios (in the K band) with group mass also 
levels off at this same point. 

The changes in M* and a are also understandable phys- 
ically in the context of group scale mergers. If groups hierar- 
chically merged but there was no change in either a or M* 
then this would be a sign that the group merger had not 
prompted galaxy scale effects. 



3.3 Multiplicity 

For the multiplicity ordered data, we will again use 10 bins 
with as near to equal numbers of galaxies as possible, but a 
3 bin version is also used to aid later analysis of results. For 
these subsets the sample is split into groups of fixed multi- 
plicity ranges: 5 < N < 9, 10 < N < 19 and N > 20. These 
3 bins contain 6516, 4393 and 6896 galaxies respectively. 

The results can be found in figure[S] In terms of M* , we 
see little systematic change for the overall population until 
we reach the very lowest multiplicities. Indeed, splitting by 
colour, the blue galaxies only have a fainter M* in the lowest 
bin. Generally, though, the trends are in line with our other 
group orderings. 

Interestingly, the lowest multiplicity systems, in this 
case all groups with only 5 members, also experience more 
dwarf suppression for blue galaxies than any of the other 
methods for sorting the data, with a ~ —0.8. This implies 
that a system being numerically poor has more of an effect 
on the lack of late-type dwarfs than the mass of the system. 
The red galaxies have the same extreme value of a ~ —0.3 
as the small groups in the other orderings. At the other end, 
the a parameters settle to their normal rich system values, 
which is understandable since multiplicity correlates very 
tightly with virial mass in the largest systems (central plot 
of figure [3J) . 

In the lower panels of figure [5] the data is split into just 
3 bins, with specific multiplicity boundaries (at 10 and 20). 
The smoothing out of values is what you might expect when 
combining composite groups. The trend for M* is still not 
prominent for these larger composites, despite the smaller 
error bars, it being slightly fainter only in the 5 < N < 
9 bin. Even the early type galaxies have a range of only 
M* ~ -19.4 to M* ~ -19.8), only just outside of the 95% 
confidence intervals, a much smaller effect than seen in mass 
ordered groups. 

The trends in a, on the other hand, are extremely sim- 
ilar to those seen for the other group orderings. The dif- 
ferences in the best fit a values for early-type galaxies are 
again quite pronounced, and are clearly seen in the lower 
panel with only three multiplicity ranges. 

The suppression of the faint-end in the low-multiplicity 
systems can also be seen by plotting a variation on the 
dwarf-to-giant ratio (see Figure [BJ, or directly from the lu- 
minosity distributions themselves (right hand panel) . In the 
former, the fraction of galaxies brighter than the limiting ab- 
solute magnitude at the high redshift extreme of the survey 
volume (Mb- = —18.25 at z — 0.1) in each group is plotted 
as a function of group multiplicity and redshift. The lowest 
multiplicity systems maintain a larger fraction of galaxies 
with Mi,. < —18.25 to lower redshifts than their larger mul- 
tiplicity counterparts, a large number still have fractional 
giant population of ~ 1 (i.e. no dwarfs) at z = 0.07. At 
such redshifts the limiting magnitude of the survey is com- 
fortably fainter than this dividing line; assuming a typical 
galaxy k-correction it will be ~ —17. 25. Thus if th e stan - 
dard 2dFGRS Schechter LF values of IColless et~aH (|200lD 
were used, the expected population fraction of giants would 
be only 0.58 at this redshift. The point at which this fraction 
consistently reaches 0.5 is only at z ~ 0.055 (i.e. when the 
limiting magnitude is —16.75) for N = 5 groups. In com- 
parison TV = 10 groups achieve this level at z ~ 0.07, as 
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Figure 4. Luminosity functions for groups from the 2PIGG catalogue as found via Maximum Likelihood for blue, red and all galaxies 
(blue, red and green colours respectively). Plots show best fit M* and a for groups sorted into 10 bins via virial mass estimates, x-axis 
errors are the standard deviation of the data, and the y-axis errors are the 95% confidence levels computed from the log-likelihood. 



expected. The overall luminosity distributions in the right 
hand plot demonstrate visually the clear differences in faint 
end slopes between small and large multiplicity groups for 
red galaxies in particular. 

3.4 Field LF and Lower Multiplicities 

Recall that the low mass groups are a more specific envi- 
ronment than the general field (where local densities are 
computed over much larger volumes than the typical group 
sizes). It is therefore not necessary (and indeed it is not the 
case) that the faint end of the LF in such groups should tend 
to that of the field. Field LFs tend to have a ~ —1.1 rather 
than the shallower slopes for small groups seen here. 

Despite the non-trivial nature of the 'field', it is instruc- 
tive to compare the trends discussed above with the LF pa- 
rameters for all objects that are not in 2PIGG groups of 
multiplicity 5 and greater but meet the volume limited cri- 
teria of our sample. For the un-grouped sample (i.e. TV = 1) 
we find M* = -19.16 ± 0.04 and a = -0.83 ± 0.02 and for 
the groups too small to meet our sample (2 < N < 4) we 
find M* = -19.34±0.06 and a = -0.71±0.04. It should be 
highlighted that the fit for N — 1 is of much better quality 
than for 2 < N < 4. Comparing these numbers to the green 
points in Figure[5]indicates the trend in a, and the more sub- 
tle trend in M* continues for lower mass groups that don't 
make our sample. The empirical LFs, along with the best 
fitting Schechter functions, can be found on the left panel of 
Figure for the data cuts described above. As a reference 
the composite LF for all group with N > 5 is shown, for 
which we find M* = -19.79 ± 0.03 and a = -1.01 ± 0.02. 
Noticeable in this plot is the excess of bright galaxies for 
the N > 5 population, despite the smaller integrated pop- 
ulation. This can be interpreted as cruder evidence of the 
variation in M* with the scale of environment (un-grouped 
through to grouped), however we are only breaking the pop- 
ulation into 3 large subsets. 



For the red population, selected as above, we find for 
N = 1 groups that M* = -18.95 ± 0.03 and a = -0.23 ± 
0.03, and for 2 < N < 4 groups M* = -19.18 ± 0.03 and 
a = —0.23 ± 0.03. Again, these values continue the strong 
trend in a into lower multiplicity/ mass groups. Finally, for 
the blue galaxies we find for N — 1 groups M* = —19.49 ± 
0.03 and a = -1.23 ± 0.02, and for 2 < N < 4 groups 
M* = -19.62 ± 0.04 and a = -1.17 ± 0.03. Neither of 
these values varies significantly from that lowest mass groups 
plotted in figure [5] and it is reasonable to conclude that any 
evolution of the blue galaxy luminosity function is weak or 
non-existant at low masses/ multiplicities. 



4 LF FOR YANG-2DF GROUPS 

The Yang-2dF group catalogue is based on the same sur- 
vey data as the 2PIGG catalogue. However, as discussed in 
Section 2, the method for grouping the galaxies is markedly 
different. It is therefore instructive to see how different the 
best fit luminosity functions are for the two grouping meth- 
ods. 

In so far as is possible the same composite grouping 
criteria have been used: a requirement of at least 5 group 
members, the same survey volume and a division based on 
the total mass of the groups - though this is now based on 
the method described by Yang et al. (2005b). Due to the 
poorer number statistics for the Yang-2dF groups the data 
has been split into 5 composite groups for each method of 
ordering, as opposed to the 10 used for 2PIGG. 

4.1 Similar Versus Non-Similar Groups 

As mentioned previously, a certain fraction of groups in the 
Yang-2dF and 2PIGG catalogues contain a large fraction of 
overlapping galaxy populations. The first obvious compari- 
son to make, therefore, is between those groups that appear 
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Figure 5. Luminosity functions for groups from the 2PIGG catalogue as found via Maximum Likelihood for blue, red and all galaxies 
(blue, red and green colours respectively). Plots show best fit M* and a. Both top and bottom plots are sorted by group multiplicity. 
The top plot has 10 splits based on maintaining galaxy count in bins. The bottom plot has three bins of roughly similar galaxy count, 
but specifically the 3 ranges: 5 < N < 9, 10 < N < 19 and N > 20. x-axis errors are the 25% and 75% quantiles, and the y-axis errors 
are the 95% confidence levels computed from the log-likelihood. 



in both catalogues and that have the closest matching group 
centre. This latter requirement is so that we can ensure the 
comparison made is between the same central haloes, and 
means in the 47 cases where multiple Yang-2dF groups sig- 
nificantly match to a 2PIGG, only the closest groups in co- 
moving space are compared. Obviously there is ambiguity 
in how best to match and compare such catalogues, but the 
aim here is to compare how the extra non-core galaxies af- 
fect the LF of the whole group- knowing that the Yang-2dF 
groups will be describing the core population only. 

The best fit LF for matched groups based on the 2PIGG 
catalogue is M* = -20.07 ± 0.05 and a = -1.14 ± 0.03, 
compared to M* = -20.28 ± 0.07 and a = -1.16 ± 0.04 for 
the matched Yang-2dF groups. The shapes of the LFs for the 
similar groups in the two catalogues agree best at the bright 



end and diverge slightly as a function of faintness, though 
evidently without a huge effect on the LF parameters as a 
large fraction of the galaxies in these two subsets will be the 
same. The biggest difference is simply in the normalisation, 
and suggests that where the two algorithms agree the LF 
is not varying greatly between the group core (reflected by 
Yang-2dF groups) and the extended infall region (captured 
by 2PIGG). 

The remaining unmatched groups have a M* = 
-19.63 ± 0.03 and a = -0.92 ± 0.02 for 2PIGG, and 
M* = -20.42 ± 0.12 and a = -1.28 ± 0.06 for Yang-2dF 
groups. It is worth noting that the best fit LFs for these sub- 
sets generally have a low p- value (< 0.1), so it is dangerous 
to consider the best fit parameters alone. The best fit LFs 
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Figure 6. Left plot shows the median fraction of galaxies with Mi,. < —18.25 in groups as a function of original multiplicity and 
redshift. It is quite clear that the lower multiplicity systems have a larger number of 100% brighter than giant limit galaxies at lower 
redshifts, indicating a lack of faint galaxies even where they would be expected to be observable in low multiplicity systems. Error bars 
are calculated from the error in the mean. Right plot show the counts-corrected luminosity function for N < 9 and TV > 20 groups. Blue 
depicts blue galaxies (as defined by the b — r cut), red depicts red galaxies. Clearly the biggest discrepancy occurs at the faint end for 
early-types, the low multiplicity systems being massively depleted. 




Figure 7. Left panel compares the LF for all 2PIGGs used in our sample, 2PIGGs that are grouped but are not in our sample by virtue 
of multiplicity (2 < N < 4) and all un-grouped galaxies. Comparison of LFs formed by groups that are matched in the 2PIGG and Yang 
group catalogues, and groups which are unmatched. Data and best fit LFs are as depicted in the legend. 



along with the binned version of the data can be found in 
Figure [7] 

A point of interest is that the matched and unmatched 
LFs for the same catalogue generally agree (or more specif- 
ically the LFs converge) at the bright and faint ends of the 
LF. In between these extremes the LFs diverge consider- 
ably. In the case of 2PIGG, the unmatched groups have 
a considerable excess of galaxies in the moderate luminos- 
ity ( — 20 < Mbj < —19) regime. In contrast, this is the 



region where the unmatched Yang-2dF groups display a 
relative deficit. The brighter M* value for the unmatched 
Yang groups is understandable since these groups (due to 
the method used to group them) will typically contain a 
dominant BCG. Since 2PIGG has no such bias, it is to be 
expected that in situations where the two catalogues fail 
to match a group the Yang groups will possess a brighter 
central galaxy. 

The unmatched group LFs are important because they 
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give an insight into the most divergent behaviour between 
the grouping algorithms. The regions that have a relative 
excess of moderate luminosity galaxies will generally be 
preferentially grouped by 2PIGG, and the reverse appears 
to be true for the Yang-2dF catalogue. This might indi- 
cate that the linking length used in 2PIGG for clustered 
moderate-luminosity galaxies is relatively generous. By this 
we mean that the global LF used for determining grouping 
over-estimates the rarity of moderate luminosity galaxies in 
dense environments, these galaxies are apparently the key 
objects that determine successful grouping. In contrast, the 
lack of this moderate luminosity excess in the non-matched 
Yang-2dF groups implies that these galaxies do not add as 
much mass to the system (as calculated by the CLF), and do 
not increase the group-scale linking length as significantly. 
Ideally no one type of galaxy should dominate proceedings 
in this manner, but it does indicate that the groups with 
dominant central galaxies or very dense regions (i.e. regions 
the Yang et al. algorithm is designed to find) are relatively 
lacking in moderate luminosity galaxies, the expected signa- 
ture of 'cannibalism'. Obviously this indicates environmen- 
tal dependence, but more generally this demonstrates that 
the Yang et al. algorithm is more prescriptive (compared to 
2PIGG) about how the groups are built. This is reflected by 
the smaller variation in LF parameters between the similar 
and non-similar Yang-2dF/2PIGG systems. 

4.2 Mass Ordered Yang-2dF Groups 

The luminosity functions ordered by total mass can be found 
in Figure [8] The main variation compared to 2PIGG would 
seem to be steeper values for a and brighter values for M* . 
The trends for M* appear to be much better defined than for 
the mass ordering of 2PIGG. However this is largely due to 
differences in how the total masses of the groups are calcu- 
lated. The 2PIGG method assumes a degree of virialisation 
and is based on the velocity dispersion and a characteristic 
radius, while the Yang et al. method is based on the lu- 
minosity of the groups and an assumed mass-to-light ratio. 
As such it is largely the same as ranking by the total gi- 
ant magnitude, as was done for 2PIGG in Robotham et al. 
(2006). 

The combined, blue and red galaxies are all significantly 
brighter in M* for all masses compared to 2PIGG (though 
with quite large error bars). The red galaxies are found to 
have much steeper best fits throughout the range of com- 
posite groups compared to 2PIGG: the range covered is 
-0.8 > a > -1.2 as opposed to -0.3 > a > -1 for 2PIGG. 
The late types show no significant variations in a, but are 
consistent with 2PIGG results given the large errors. The 'all 
galaxies' LFs show a modest steepening at all masses, sug- 
gesting there is notable variation in how early-type galaxies 
are treated. 



4.3 Multiplicity Ordered Yang-2dF Groups 

As for 2PIGG, the other method used for ordering the data 
was by group multiplicity. The results for this can be found 
in Figure [5] the top plots show the best fits for 5 compos- 
ite groups with efforts to maintain galaxy numbers in each 
composite group, the bottom plots use the same fixed mul- 
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Table 1 . Table comparing the numbers of different galaxy types 
(red/blue) for 2PIGG and Yang-2dF groups for different multi- 
plicity subsets. 

tiplicity limits as previously used for 2PIGG (5 < N < 9, 
10 < N < 19 and N > 20). 

For the 5 bin split via multiplicity, the errors for the 
late- types are typically very large, so the apparent faintening 
trend as a function of multiplicity may not be reliable. The 
change in M* for early-types for the lowest multiplicities 
again appears to be present in this data, though as with all 
the Yang-2dF composites, it takes a brighter value than for 
2PIGG. 

The best fit values for a are again generally steeper 
than for 2PIGG for equivalent composite group subsets, but 
the trends are much the same. The early-type galaxies show 
very clear steepening as a function of multiplicity and, much 
like was found for 2PIGG, the low multiplicity systems show 
the most dramatically depressed faint end with a ~ —0.7. 
The blue galaxies have a ~ —1.5 throughout. 

4.4 LF Differences Between Grouping Algorithms 

Whilst there is a good deal of agreement between 2PIGG 
and Yang-2dF group luminosity functions, it is the differ- 
ences that are probably the more interesting. As discussed 
above, the most distinct variation appears to be that a val- 
ues are routinely steeper for the Yang-2dF groups than for 
2PIGG. This is worth investigating in more detail. 

Figure [TO] is a like for like comparison between the two 
catalogues at the small (5 < N < 9) and large (20 < N) 
extreme for red galaxies, and it is immediately clear that 
whilst the two populations display differences for both of 
these multiplicity subsets, it is in the lowest multiplicity cut 
that the most significant discrepancies appear. 

Considering the large multiplicities first, the Yang-2dF 
group data is more similar to 2PIGG in its LF than the 
best fit implies, as the underlying LF is not well reproduced 
by the fitted Schechter function (as reflected by a very low 
p- value for the best fit of 0.0014). The major disagreement 
between the two LFs is at the bright end, where the number 
of giants in the Yang groups is comparitively high and thus 
the dwarf-to-giant ratio is lower. Elsewhere the results indi- 
cate a similar PDF for the LF but a less rich population - for 
most bins there is a factor ~ 3 difference. From the variation 
in gradient seen, it appears that the Yang-2dF groups would 
be better fit by a combination of separate giant and dwarf 
(satellite) populations (c.f. Popesso et al. 2005, for clusters). 
Specifically, the gradient steepens at M&. ~ —20.5, when for 
a pure Schechter function it will always decrease. 

The 5 < N < 9 LFs show a great deal of disagreement 
between —20 < Mb < —18, where 2PIGG displays a general 
excess of moderate luminosity galaxies compared to Yang- 
2dF groups. Interestingly, the two populations converge al- 
most completely at the giant and faint extremes, suggesting 
that for these low multiplicity systems the most bound struc- 



The Luminosity Function in Galaxy Groups 13 



-*- 



13 



T" 



I 

14 



T" 



13,5 14 14.5 

log 10 (Group Mass / M su „) 



+ 



13 



T 



14 



T" 



13,5 14 14.5 

log 10 (Group Mass / M su „) 



Figure 8. Luminosity functions for Yang-2dF groups as found via Maximum Likelihood for blue, red and all galaxies (blue, red and 
green colours respectively). Plots show best fit M* and a for groups sorted into 5 bins via total mass estimates, x-axis errors are the 
standard deviation of the data, and the y-axis errors are the 95% confidence levels computed from the log-likelihood. 



ture is at either extreme of the luminosity function, and the 
tenuous structure (as grouped by 2PIGG) tends to be re- 
flected by the moderate luminosity red galaxies. Physically, 
2PIGG may be 'attaching' extra moderate luminosity ob- 
jects which are/were the central galaxy of their own halo 
when grouping objects, whereas the Yang et al. algorithm 
does not. 

Figure [11] displays the different LFs for blue galaxies 
using the same multiplicity subsets. The blue galaxies in 
Yang-2dF groups suffer poor number statistics for both mul- 
tiplicity cuts, as implied by the larger error bars, and subse- 
quently have a noisier distribution around the best fit. From 
Table 1, the lower multiplicity Yang-2dF groups are substan- 
tially less rich in blue galaxies than their 2PIGG counter- 
parts (34% of the galaxies compared to 53%) and the high 
multiplicity systems are relatively even poorer in blue galax- 
ies (19% compared to 38%). Thus the ratio of red to blue 
galaxies is always larger for Yang-2dF groups, but this ratio 
increases dramatically as a function of multiplicity, indicat- 
ing that in more bound systems (in this case, the Yang-2dF 
groups) a greater fraction of blue late-types (dominating the 
blue sample) have made the transition to early-types (dom- 
inating the red sample), even at fixed multiplicity. As these 
more bound g roups wil l have higher density, this then re- 
flects the usual Dressier ( 1980) morphology-density relation, 
but within fairly small groups. 

In addition, the blue galaxy LF remains fairly consis- 
tent in its PDF (i.e. shape), between the catalogues, whilst 
the red galaxy LF is altered as a function of both richness 
and the bounding criterion used. The halo-based method 
for the Yang-2dF groups clearly produces a relative deficit 
in moderate luminosity red galaxies compared to the FoF 
2PIGG, particularly for low multiplicity systems. This may 
indicate that early-type galaxies possess a LF shape that 
varies significantly as a function of environment, i.e. environ- 
ment affects different luminosity galaxies quite differently, 



whilst late-type LFs mostly vary in richness as a function of 
environment, the shape remaining relatively fixed. 

There is a caveat to making strong conclusions regard- 
ing the discrepancy between 2PIGG and Yang-2dF groups 
early /late- type ratios, since it can be partially explained by 
a relative excess of interlopers in 2PIGG systems. These in- 
terlopers are likely to be field galaxies, and will tend to be 
late-type (blue) galaxies. Whilst this might partially explain 
the discrepancies in relative richness, it does not naturally 
explain variations in LF shape seen for the red galaxy sam- 
ple. 



5 POSITION DEPENDENCE 

To explore this a little more we consider the variation of 
the luminosity function with position within 'halos' for the 
2PIGG catalogue. Taking the centres determined previously, 
we have spatially 'stacked' the groups in the three broader 
multiplicity composites (to maximise the galaxy numbers) 
and obtained LFs for the combined and the separated blue 
and red populations, as defined above. 

The LF was calculated for both types within a projected 
radius of 0.5 r rm3 , and outside of 1 r rrns for each multiplicity. 
These limits were chosen so as to make sure a large fraction 
of galaxies in these areas are physically within the 3D radii 
regimes of interest. The results for these fits can be found 
in tabled 

We can see that the outer regions maintain a similar 
LF shape through all environments and galaxy types, within 
95% confidence errors. The central LF for late- types does ap- 
pear to change slightly as a function of multiplicity, however 
this variation is only just outside the confidence range for the 
parameters. The most dramatic variation is certainly that 
seen for early-type (red) centrally located galaxies. The key 
parameter of interest, since it will determine the majority 
of galaxy counts in a given volume, is a for the early-type 
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Figure 9. Luminosity functions for Yang-2dF groups as found via Maximum Likelihood for blue, red and all galaxies (blue, red and 
green colours respectively). Plots show best fit M* and a. Both top and bottom plots are sorted by group multiplicity. The top plot has 
5 splits based on maintaining galaxy count in bins. The bottom plot has three bins of roughly similar galaxy count, but specifically the 3 
ranges: 5 < TV < 9, 10 < N < 19 and N > 20. x-axis errors are the 25% and 75% quantiles, and the y-axis errors are the 95% confidence 
levels computed from the log-likelihood. 



galaxies. For the multiplicity 5-9 groups it changes from 
a = —0.52 centrally to a = —0.21 in the outer parts (the 
latter remaining similar for the other multiplicities), but the 
central LFs have a = —0.85 for groups with 10-19 members 
and a — —1.00 for the largest systems. In the latter case the 
LF is actually quite a poor fit (p- value of 0.01), and a visual 
inspection of the raw LF indicates a two component model 
is almost certainly more appropriate, but regardless of this, 
the faint end is clearly markedly steeper for the 10-19 and 
20+ member groups compared to 5-9 multiplicities. 

Following ICroton et al.l (|2005l ). this implies that the 
outer regions of groups are very similar environments and 
apparently influenced almost uniquely by the local (instan- 
taneous) density of the environment, which will be similar at 
the outer limits of all FOF grouped system (see also the dis- 



cussion in. Baldrv at al.| [2006). This has the added implica- 
tion that beyond r rms , interactions with the denser cores of 
groups have had little impact on the galaxy population, else 
we would expect to see some variation in systems with much 
denser central regions. These results are consistent with the 
work of, for example, Haines et al. (2006), who found that 
the dwarf population demonstrated a sharp transition be- 
tween early and late-type dominance near 1 virial radius 
from the centre of Abell 2199. The conclusion of that work 
was that dwarf galaxies only become passive once they have 
become satellites within a more massive halo. More gener- 
ally, Popesso et al. (2006) suggest that galaxy transforma- 
tion occurs at around 0.7r2oo, and it would appear that sim- 
ilar conclusions can be drawn on scales much smaller than 
the clusters used in their work. Such a picture is consistent 
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Figure 10. Red galaxy LFs for Yang-2dF groups and 2PIGG. Shown are the completeness corrected counts and the best fit maximum 
likelihood LF. The low multiplicity (5 < N < 9) comparison are shown in the left plots, the large multiplicity (20 > N) comparisons 
on the right. Disagreement between the data and best fit is partially caused by the use of maximum likelihood (not binned data) and 
because the mid-point of the bin is used for scaling purposes. 
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Figure 11. Blue galaxy LFs for Yang-2dF groups and 2PIGG. Shown are the completeness corrected counts and the best fit maximum 
likelihood LF. The low multiplicity (5 < N < 9) comparison are shown in the left plots, the large multiplicity (20 > N) comparisons 
on the right. Disagreement between the data and best fit is partially caused by the use of maximum likelihood (not binned data) and 
because the mid-point of the bin is used for scaling purposes. 



wit h semi-analytic models of star formation suppression (see 
e.g. iBower et al]l200d : ICattaneo et"ai1l2007l ). 

It is important to note that r rms appears to be the 
important threshold for defining the region where the LF 
dramatically changes shape, most significantly for early-type 
galaxies. Outside of this, grouped objects appear to share 
a largely similar LF, regardless of the size of the system 
(as determined by multiplicity). The strength of this effect 
is readily observable in Figure 1121 which compares a as a 
function of the scaled radius and the actual projected radius. 



It seems clear that the scaled radius is a better tool for 
finding regions with similar LFs, and all three environments 
used in this work are most similar beyond 1.0 r rmB . 



5.1 Luminosity and Distance Ranks 

Related to these positional variations is the question of 
luminosity-distance correlations. If a group is well virialised 
and evolved, then the process of merging will tend to place 
the most luminous galaxies closest to the central galaxy. 



16 A. Robotham et al. 




o 
d 



d 



CO 

d 
i 



7 









J- <><l> J . 










i 



Projected Radius/ Scale Radii 



Projected Radius/ Mpc 
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projected radius r. In all cases the divisions are based on maintaining galaxy counts in the subset (hence similar 95 confidence regions), 
and the midpoints of the bin is used for plotting purposes. 



The definition of the group centre varies between 2PIGG 
and Yang-2dF groups. 2PIGG uses an iterative process to 
continually remove the galaxy furthest from the complete- 
ness weighted centre until only two are left. The final step 
is to select either the galaxy with the largest weight, or if 
identical the galaxy with the most bj flux. The Yang al- 
gorithm uses a luminosity weighted group centre from the 
outset, so at least the brightest member will have an en- 
hanced tendency to be near to the defined group centre. We 
are interested in uncovering the luminosity gradient from 
the brightest group member rather than relative to a some- 
what arbitrarily defined group centre. With this in mind all, 
group centres were redefined from the brightest galaxy for 
both 2PIGG and Yang groups, and biases caused by differing 
definitions of the group centre were therefore removed. 

As a simple measurement, a real merger gradient im- 
plies a negative correlation between galaxy luminosity and 
distance from the brightest galaxy in the group (or a posi- 
tive correlation if we use magnitudes) . To remove the issue of 
group scaling, the Pearson product-moment coefficient was 
calculated for each group in turn (from both the 2PIGG 
and Yang-2dF catalogues), and this was compared to the 
distribution expected for randomly generated data with the 
same multiplicities. The form of the Pearson statistic (P r ) 
measured was 



1 



/ M rt - M r 



= 1 ^ 



(9) 



where M rp is the absolute tf band magnitude and r is the 
projected distance from the brightest galaxy for each galaxy. 
This statistic returns if the correlation is non-existant, 1 if 
it is perfectly linear and positive and —1 if it is perfectly lin- 
ear and negative. All distances are relative to the brightest 
group member so the brightest galaxy is removed from the 
calculation, the remaining gradient is for the second rank 



galaxy and fainter. Since it is the projected distance that 
is used the correlation, if present, will be reduced in sig- 
nificance by the random positions of group foreground and 
background galaxies. 

The results for the Pearson correlation can be found in 
the left hand plot of Figure 1131 which describes the PDFs 
for the P r statistic for both 2PIGG and Yang-2dF data com- 
pared to the distribution expected for a randomly generated 
sample of equal multiplicity groups with randomised magni- 
tudes and distances. It is immediately clear from these plots 
that whilst the 2PIGG data shows no evidence for any cor- 
relation, the Yang-2dF group data shows a strong tendency 
towards positive correlation. 

This difference presumably arises from the way that the 
groups were originally constructed. In the case of the 2PIGG 
FOF algorithm there is no bias towards more massive galax- 
ies being grouped together, the linking treats all galaxies 
equally at a given redshift. The Yang algorithm begins with 
progenitor groups constructed from highly overdense regions 
in space or luminous single galaxies, and such a prior will 
clearly generate groups with dominant central galaxies. The 
apparent presence of positive correlation between distance to 
the central galaxy and the magnitude indicates that Yang- 
2dF groups show evidence for some sort of merger history, 
and thus describe a more evolved and virialised structure 
than 2PIGGs. This is consistent with the results from pre- 
vious sections which also imply the Yang-2dF groups to be 
more tightly bound and virialised systems. 

As a point of comparison, exactly the same methodol- 
ogy was applied to groups found in the Millenn ium Simula- 
tion (|Springel et al.ll2005l : be Lucia etaH 120061 ) . Again the 
correlation was calculated between the projected distance 
from the central galaxy and the galaxy r-band magnitude, 
where the projection was chosen randomly. The results for 
this simulated data can be found in the right-hand plot of 
Figure [131 Immediately obvious is the extreme bias towards 
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multiplicities are sampled according to the group occupation. Right-hand plot is the same but generated for a sample of Millennium 
Simulation groups. 



the positive correlation gradients expected for evolved and 
bound systems, with a significance even higher than for the 
Yang groups. KS-test statistics were calculated between all 
real and mock P r distributions and can be found in the rel- 
evant plots of Figure [Tol Whilst the Yang-2dF groups reject 
a purely random distribution of galaxy magnitudes relative 
to the central galaxy at the 99% level, the Millennium Sim- 
ulation groups reject this at above the 7a level. The 2PIGG 
sample demonstrates an effectively random distribution of 
magnitudes within the groups and has a KS-test p-value 
of 0.44. This adds further weight to the interpretation that 
Yang-2dF groups describe more evolved systems, similar to 
those generated using semi-analytic modelling in the Mil- 
lennium Simulation. De Lucia et al. (2006) interpret the 
simulation results as implying that the radial distance cor- 
relates with the time at whic h galaxies were accreted onto 
the central structure (see also lGao et al.lr2004l ). 



6 SUMMARY 

In this paper a detailed comparison between the grouping 
characteristics of the 2PIGG and Yang group catalogues 
was made, with particular reference to the group luminosity 
functions and how these are affected by the different group- 
ing algorithms. We have extended the work of Eke et al. 
(2004a) and Robotham et al. (2006) on the luminosity func- 
tion of galaxies in groups of different masses by (a) using a 
maximum likelihood method for determining the LF param- 
eters, (b) using different indicators of the group halo mass, 
viz. a mass direct virial mass estimate and the number of 
galaxies in the group. 

For 2PIGG, virial masses were calculated as in Eke et al. 
(2004b) from the velocity dispersion and the rms distance 
from the group centre of the member galaxies . The same 
trends as previously found in iRobotham et all |2006l ) are 



still generally present, but are somewhat less well defined 
for M* , especially in the case of blue galaxies. The faint 
end slope, a, for the total population changes in a similar 
way to that seen in previous work, but the subsets behave 
significantly differently: red galaxies now showing a large 
steepening, blue ones little or no change until we reach the 
highest masses sampled. 

In terms of ordering by multiplicity, it appears that only 
the groups with very few (giant) members differ substan- 
tially in M* from the larger objects. This is also the case 
for a, though the blue galaxies in low multiplicity groups 
have the least negative slope of any blue galaxies sampled. 
This appears to imply that a system being numerically poor 
has more effect on the (lack of) blue dwarfs than the group's 
mass. This effect is also seen if we consider the relative num- 
bers of galaxies brighter than and fainter than Mt = —18.25 
(for groups sufficiently nearby for the latter to be visible). A 
surprising number of the lowest multiplicity groups have few 
or no galaxies in the fainter range (typically Mi, between 
-18.25 and -17.25). 

The Yang-2dF groups provide a useful counterpoint to 
the 2PIGG groups as they were constructed via an algo- 
rithm which essentially works via the binding of a galaxy to 
a pre-existing group core. In general terms, the Yang-2dF 
groups demonstrate the same LF parameter trends as be- 
fore. However, there are a number of distinct differences in 
detail. 

Firstly, even for groups clearly centred at the same posi- 
tion in the two catalogues (i.e. representing the same halo), 
the LF shapes appear to diverge somewhat for the fainter 
galaxies, though the LF parameters are closely similar. Sec- 
ond, the 2PIGG groups not closely matched by a Yang-2dF 
group have a fainter M* and flatter a than the rest, whereas 
unmatched Yang-2dF groups have brighter M* and steeper 
a, though we should also note that the LFs for these subsets 



18 A. Robotham et al. 



are not, in fact, well fitted by a Schechter function. A two 
component fit would be preferred for the Yang-2dF groups. 

In terms of the LF shapes, it appears that the matched 
and unmatched group LFs differ primarily through an ex- 
cess (shortfall) of intermediate luminosity galaxies in the 
unmatched 2PIGG (Yang-2dF) groups. These unmatched 
groups are important as they are the ones which most clearly 
reflect the behaviour of the grouping algorithms. It appears 
that 2PIGG will quite generously link objects where there 
are numerous medium luminosity galaxies, whereas the Yang 
algorithm docs not (presumably because they do not add 
sufficient mass to make the group encompass further ob- 
jects). 

Physically, this indicates that groups with dominant 
central galaxies or dense central cores, of the sort the Yang 
algorithm is designed to find, are relatively deficient in 
medium mass objects, consistent with them being more dy- 
namically evolved. 2PIGG, on the other hand, may tend 
to link subgroups, adding additional moderate mass objects 
which are/were the central galaxy of their own (infalling) 
halo. 

Among the Yang-2dF groups themselves, while the er- 
rors are larger than for 2PIGG, the mass (here really just 
a scale version of luminosity) ordered composites show the 
same trends as the 2PIGG groups, with M* brightening and 
a steepening for the more luminous groups. However, as 
implied by the comparison above, the actual values differ, 
sometimes substantially. In particularly M* is ~ 0.5 mag- 
nitudes or more brighter at all masses for Yang-2dF groups 
than 2PIGG groups. Red galaxies have consistently steeper 
a in the Yang-2dF groups, while, given the errors the blue 
galaxies are consistent with the same slopes as those in 
2PIGG (or, indeed, with no change with group mass). 

The multiplicity ordered Yang-2dF groups show the 
same changes relative to 2PIGG multiplicity ordered groups 
as seen for the mass ordered groups, and as in 2PIGG the 
least negative a is seen for red galaxies in low multiplicity 
groups, blue galaxies having a similar a at all multiplicities. 
Thus the Yang-2dF groups contain a larger fraction of red 
early type galaxies at the same group mass or multiplicity 
than seen in 2PIGG and this difference is largest for the 
highest multiplicity groups. Generally, this deficit of low- 
mass passive galaxies can be directly interpret as a build up 
of the red sequence as a function of halo mass, and indirectly 
as a build up of the red sequence with redshift. 

LFs calculated for galaxies within 0.5 r rma and outside 
1 r rms show that the outer regions have similar LF shapes 
for all group multiplicities for both red and blue galaxies. 
The central LF for late type galaxies changes slightly with 
multiplicity, but the major change is for centrally located 
early type galaxies; a is significantly steeper than for the 
outer regions at all multiplicities and steepens even more for 
higher multiplicities. Looked at in a slightly different way, 
we have seen that the Yang-2dF groups tend to be denser 
and probably more evolved, and the change in population 
with radius for such systems is reflected in a strong rank 
correlation between luminosity and group-centric distance 
within these groups. 

Thus we should finally conclude that the distribution of 
luminosities of galaxies in groups depends both on the mass 
of the group, as determined via any suitable (and consis- 
tent) manner and on the (density) structure of the group, 



so that different grouping algorithms will generate different 
luminosity distributions for the grouped galaxies. 

Finally distance-rank magnitude relations were consid- 
ered. Only the Yang groups demonstrated any evidence of a 
correlation between a galaxy's position relative to the cen- 
tre and its luminosity (i.e. brighter galaxies have a weak 
tendency to be nearer the group centre). 2PIGG possessed 
no such gradient, the conclusion being the non-prescriptive 
nature of luminosity ordering in the FOF algorithm destroys 
any weakly observable trend. A comparison was made to the 
semi-analytic galaxy catalogue of De Lucia et al. (2006) built 
on top of the Millennium Simulation. Interestingly this the- 
oretical work demonstrated much stronger gradients than 
present even in Yang-2dF groups, and may indicate a semi- 
analytic formalism that is too stringent. 
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Environment 


M* 


OL 


p- value 


low M* 


high M* 


low oc 


high a 


5-9 Blue Cen 


-19.93 


— 1.17 


0.91 


-20.11 


-19.76 


-1.27 


-1.06 


5-9 Blue Out 


-19.71 


— 1.17 


0.44 


-19.82 


-19.61 


-1.22 


-1.08 


5-9 Red Cen 


-19.75 


-0.52 


0.69 


-19.86 


-19.64 


-0.62 


-0.43 


5-9 Red Out 


-19.10 


-0.21 


0.35 


-19.18 


-19.03 


-0.33 


-0.14 


10-19 Blue Cen 


-20.09 


-1.32 


0.84 


-20.31 


-19.89 


-1.43 


-1.21 


10-19 Blue Out 


-19.59 


-1.17 


0.78 


-19.75 


-19.41 


-1.28 


-1.04 


10-19 Red Cen 


-20.05 


-0.85 


0.61 


-20.17 


-19.94 


-0.93 


-0.77 


10-19 Red Out 


-19.02 


-0.35 


0.99 


-19.16 


-18.89 


-0.51 


-0.19 


20+ Blue Cen 


-20.19 


-1.41 


0.21 


-20.36 


-20.03 


-1.49 


-1.33 


20+ Blue Out 


-19.59 


-1.17 


0.73 


-19.81 


-19.37 


-1.31 


-1.00 


20+ Red Cen 


-19.95 


-1.00 


0.01 


-20.03 


-19.86 


-1.06 


-0.95 


20+ Red Out 


-19.08 


-0.21 


0.89 


-19.23 


-18.94 


-0.45 


-0.08 



Table 2. LFs for inner and outer regions of groups. 



